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Abstract. Quantum confined devices that manipulate single electrons in graphene are emerging 
as attractive candidates for nanoelectronics applications [1^]. Previous experiments have employed 
etched graphene nanostmctures, but edge and substrate disorder severely limit device functionality [5, 
6]. Here we present a technique that builds quantum confined structures in suspended bilayer graphene 
with tunnel barriers defined by external electric fields that break layer inversion symmetry [7,8], thereby 
eliminating both edge and substrate disorder. We report clean quantum dot formation in two regimes: 
at zero magnetic field B using the single particle energy gap induced by a perpendicular electric field 
and Bi B > using the quantum Hall ferromagnet 1^ = gap for confinement [9-11]. Coulomb 
blockade oscillations exhibit periodicity consistent with electrostatic simulations based on local top gate 
geometry, a direct demonstration of local control over the band structure of graphene. This technology 
integrates single electron transport with high device quality and access to vibrational modes, enabling 
broad applications from electromechanical sensors to quantum bits. 

Nanopattemed graphene devices, from field-effect transistors to quantum dots [2,3, 12], have been 
the subject of intensive research due to their novel electronic properties and two-dimensional struc- 
ture [13, 14]. For example, nanostructured carbon is a promising candidate for spin-based quantum 
computation [1] due to the ability to suppress hyperfine coupling to nuclear spins, a dominant source 
of spin decoherence [15-17], by using isotopically pure ^^C. Graphene is a particularly attractive 
host for lateral quantum dots since both valley and spin indices are available to encode information, a 
feature absent in GaAs [18-20]. Yet graphene lacks an intrinsic bandgap [14], which poses a serious 
challenge for the creation of such devices. Transport properties of on-substrate graphene nanostmc- 
tures defined by etching [2, 3] are severely limited by both edge disorder and charge inhomogeneities 
arising from ionized impurities in gate dielectrics [5, 6]. The absence of spin blockade in etched 
double dots is perhaps symptomatic of these obstacles [21,22]. Unzipping carbon nanotubes yields 
clean nanoribbon dots, but this approach cannot produce arbitrarily shaped nanostructures with tun- 
able constrictions [4]. Here we report fully suspended quantum dots in bilayer graphene with smooth, 
tunable tunnel barriers defined by local electrostatic gating. Our technique, which artificially modifies 
the bandgap of bilayer graphene over nanometer scales, achieves clean electron confinement isolated 
from edge disorder. 

Bernal stacked bilayer graphene is naturally suited for bandgap control because of its rich sys- 
tem of degeneracies that couple to externally applied fields. At B = 0, breaking layer inversion 
symmetry opens an energy gap tunable up to 250 meV with an external perpendicular electric field 
E [7,8, 23-26] that can be used for confinement. In devices with low disorder and at high magnetic 
fields, gapped states emerge from Coulomb-driven effects that break its eightfold degeneracy (spin, 
valley, and orbital), resulting in quantum Hall plateaus at all integer multiples of e^/h for electron 
charge e and Planck's constant h [27]. Due to the Pauli exclusion principle. Coulomb repulsions be- 
tween electrons favor spontaneous spin and/or valley polarization (or combinations of those), known 
as quantum Hall ferromagnetism, resulting in a gap at zero carrier density that far exceeds the Zeeman 
splitting energy gfisB [9, 10]. The large exchange-enhanced energy gap of A = 1.7 meV/T mea- 
sured for the = state is ideally suited for quantum confinement [11]. Because valley and layer 
indices are identical in the lowest Landau level, one may additionally induce a tunable valley gap in 
the density of states by applying a perpendicular E field that breaks layer inversion symmetry [28]. 
This coupling of valley index to E field is the key property that enables direct experimental control 
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of the relative spin and valley gap sizes in magnetic field. 

We fabricate fully suspended quantum dots with 150 to 450 nm lithographic diameters as illus- 
trated schematically in Fig. la. Following mechanical cleavage of highly ordered pyrolytic graphite 
crystals, graphene is deposited on a 300 nm thermally grown Si02 layer, which covers a doped silicon 
substrate functioning as a global back gate. After using electron-beam lithography to define Cr/Au 
electrodes, an evaporated Si02 spacer layer, and local top gates over selected bilayers, we immerse 
the samples in HF to dissolve away the dielectric on either side of the flake. This leaves both the 
graphene and the top gates suspended (Fig. la,b and Supplementary Fig. SI). Before measurement, 
the devices are current annealed in vacuum to enhance quality. The high quality of our suspended 
flakes is evident from the full lifting of the eightfold degeneracy in the quantum Hall regime and large 
resistances attained by opening the E field induced gap at 5 = and E = 90V/nm, a hundred times 
greater than reported for on- substrate bilayers at similar electric fields [23, 28]. Measurements are 
conducted in a dilution refrigerator at an electron temperature of 1 10 mK, as determined from fits to 
Coulomb blockade oscillations. 

At B = 0, the electric field effect in bilayer graphene enables the production of quantum confined 
structures with smooth, tunable tunnel barriers defined by local gating [18], thus avoiding disor- 
der arising from the physical edge of the flake. Broken layer inversion symmetry opens a bandgap 
A oc E = (aVt — f3Vt,)/2eeo, where Vt and are top and back gate voltages with coupling factors 
a and /?, respectively, and eo is vacuum permittivity. Coupling to the back gate /? is extracted from 
Landau fans in the quantum Hall regime and the relative gate coupling a//3 can be determined from 
the Dirac peak slope in a vs. V5 plot of conductance at S = 0. Properties of individual quantum 
point contacts are described in greater detail in the Supplementary Information, where pinch-off and 
behavior consistent with conductance quantization are observed (Supplementary Fig. S2 and S3). 
Quantum dot formation at S = is illustrated schematically in Fig. Ic. To create tunnel barriers 
beneath the top gates, we induce a bandgap by applying a field E while fixing Vt and at a ratio 
that maintains zero carrier density n, where n = aVt + /JV^. In the non top-gated regions, there is a 
finite charge accumulation due to an uncompensated back gate voltage. For gates in a quantum dot 
geometry, this restricts electron transport to resonant tunneling events through the constrictions. Pe- 
riodic Coulomb blockade oscillations are observed at S = which couple to both top and back gates 
(Fig. 2a). A peak in the 2D Fourier transform corresponding to an oscillation spacing of ^ llmV 
in Vb reflects this strong periodicity (Fig. 2b), and the appearance of higher harmonics reveals the 
non-sinusoidal nature of the Coulomb blockade peaks when /c^T ^ Ec, where /c^ is Boltzmann's 
constant, T is temperature, and Ec is the dot charging energy. Coulomb diamonds shown in Fig. 
2c have symmetric structure that suggests equal tunnel coupling to both the source and drain leads. 
The dot charging energy extracted from the DC bias data is Ec ^ 0.4 meV. Fig. 2d indicates that 
the periodic Coulomb blockade oscillations have comparable capacitive coupling to each pair of top 
gates. 

Coulomb blockade oscillations can also be generated at finite B field using the exchange-enhanced 
u = gap. Here the bilayer is naturally in a gapped quantum Hall state at zero density, where high 
resistances due to quantum Hall ferromagnetism make this system ideal for confinement. An iso- 
lated puddle of charge is created by fixing the Fermi energy in the top-gated regions at the middle of 
the z/ = gap while allowing occupation of higher Landau levels elsewhere, shown schematically 
in Fig. 3a. It should be noted that measurements in the quantum Hall regime are conducted in the 
valley-polarized = state, far from the transition to the spin-polarized phase [28]. Fig. 3b shows 
over forty consecutive Coulomb blockade oscillations generated at 5.2 T in a 2-gate dot with a 400 
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nm lithographic diameter. The slopes of the resonances indicate symmetric coupling to the two top 
gates, as expected for a centrally located dot. As top gate voltages are swept to more positive values, 
peak amplitude is suppressed, revealing moderate tunnel barrier tunability. Also seen in Fig. 3b are 
interruptions in the conductance resonances (vertical and horizontal features) that couple exclusively 
to a single top gate; due to their sparse and aperiodic nature, we believe that they represent charging 
events below the gates. Coulomb blockade oscillations are robust over a wide voltage range: Fig. 
3c shows an additional forty peaks generated under new gate conditions. The Coulomb diamonds 
exhibit symmetric tunnel coupling to source and drain leads and a dot charging energy of Ec ^ OA 
meV. The strongly periodic nature of the oscillations is evident in the Fourier transform of the data 
(Fig. 3d). See Supplementary Fig. S4 for additional Coulomb blockade data. 

To demonstrate geometric control over dot size, we examine the correspondence between top gate 
dimensions and Coulomb blockade peak spacing. Measurements were performed on five dots with 
lithographic diameters ranging from 150-450 nm at magnetic fields of to 7 T. The ability to de- 
crease peak spacing by increasing lithographic dot size is illustrated in Fig. 4a-c. Fig. 4a and 4b show 
Coulomb blockade peak conductance as a function of back gate voltage V5 observed in device Dl at 
B = 5.2 T, and D2 at B = 7 T, respectively (see Supplementary Information for sample labeling 
key). Black points represent data and the red lines are fits used to extract peak positions. Fig. 4c shows 
relative peak position, V{p) — Vq, plotted as a function of peak number ]9 for the first 10 peaks of Fig. 
4a and 4b, where V{p) is the position of peak p in back gate voltage, and Vq is the position of the first 
peak. Each data set is accompanied by a corresponding plot of y{p) = [| J2q=o + 1) ~ V{q)]p^ 
where p and q are peak index numbers, representing the average peak spacing (black lines in Fig. 
4c). The dot area extracted from quantized charge tunneling is given by A = l/(/3 • AV5), where 
AV5 is the back gate voltage needed to increase dot occupation by one electron (Supplementary Ta- 
ble 1). A comparison of measured dot diameter, d = ^^/AJt^, with effective lithographic diameter, 
diii]^ = 2^ Aiith/n, indicates that d generally exceeds duth- This is contrary to the reduced dimen- 
sions observed in GaAs dots, where smaller dimensions are observed due to depletion [29]. 

To obtain a better quantitative understanding of the discrepancy between lithographic and mea- 
sured dot sizes, we use a commercial finite element analysis simulation tool (COMSOL) to calculate 
the expected dot area for each top gate geometry. The spatial carrier density profile is modeled for a 
fixed top gate potential by solving the Poisson equation assuming a metallic flake in free space (Fig. 
4d and Supplementary Fig. S5). This assumption is justified by local compressibility measurements 
of the z/ = state yielding dfi/dn = 2 x 10~^^eVni? at 2 T, which translates to a screening of 99% 
of the applied V5 voltage by the bilayer [11]. Surprisingly, one may calculate dot size purely from 
gate geometry without relying on measured gap parameters. Assuming that charge accumulation in 
the quantum dot occurs when the carrier density exceeds a fixed cutoff do, the dot size is defined as 
the area bounded by the intersection with the density distribution f{x^y) with the cutoff. Assuming 
that the tunneling probability into the dot decays exponentially with barrier width, placement of the 
cutoff at the saddle points of the density profile within the constrictions enables maximal tunneling 
without loss of confinement. The simulated dot area from this method, plotted in Fig. 4c (inset), is 
simply the area bounded by the closed contour of /(x, y) at fixed density d^ (Supplementary Fig. 
S5). Alternatively, one may calculate dot size by imposing a cutoff equaling the measured gap width 
(Supplementary Table 2) and accounting for density offsets due to a displacement of the measurement 
voltage from the charge neutrality point (see Supplementary Information). Remarkably, the cutoffs 
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extracted by these two models are equivalent to within Sn ^ lO^^cm ^, the carrier density fluctua- 
tions due to disorder in our suspended bilayers (Supplementary Fig. S5) [9]. 

Our model establishes a quantitative link between measured dot size and lithographic geometry and 
therefore may serve as a design tool for future bilayer nanodevices requiring submicron spatial con- 
trol. These include double dot systems which form the basis of a spin-based quantum computer [1]. 
The production of suspended graphene quantum dots also enables study of coupling between quan- 
tized electronic and vibrational degrees of freedom [30,31], with potential applications to nanoelec- 
tromechanical devices and the detection of quantized mechanical motion in a membrane [32-35]. 
Furthermore, the combination of high sample quality with local gating enables study of edge modes 
that emerge at the interface of broken symmetry quantum Hall states in an environment well-isolated 
from edge disorder. 

Methods 

Following mechanical exfoliation of highly oriented pyrolytic graphite crystals, graphene is de- 
posited on a 300 nm thermally grown Si02 layer, which covers a doped silicon substrate functioning 
as a global back gate. Bilayer flakes are identified based on contrast to the substrate with an optical 
microscope and later verified through quantum Hall data. Cr/Au (3/100 nm) electrodes are defined 
on selected bilayers using electron beam (ebeam) lithography, thermal evaporation, and liftoff in ace- 
tone. A Si02 spacer layer approximately 150 nm thick is deposited with ebeam evaporation after a 
second lithography step. Local top gates are placed over the Si02 spacers in a two step ebeam lithog- 
raphy process. First small features that define the tunnel barriers and constrictions are patterned using 
Cr/Au of thickness 3/75 nm, and thicker support structures constructed of 3/300 nm of Cr/Au that 
traverse the evaporated Si02 step are deposited immediately afterwards. The devices are immersed 
in 5:1 buffered oxide etch for 90 s and dried in a critical point dryer, which leaves both the graphene 
and the top gates suspended. 
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Figure legends 

Figure 1| (a) Schematic cross-section of a suspended gate-defined bilayer graphene quantum 
dot. Graphene is deposited on a 300 nm thermally grown Si02 layer, followed by electron-beam 
lithography steps to define Cr/Au electrodes, an evaporated Si02 spacer layer, and local top gates. 
Etching the dielectric on either side of the flake leaves both the flake and the top gates suspended. 
The electric field and carrier density profiles are controlled with back and top gate voltages V5 and Vt, 
while application of a bias Vgd across the electrodes enables transport measurements, (b) Scanning 
electron micrograph of quantum dot device similar to D4. Bilayer graphene (not visible) is suspended 
between two electrodes below local top gates. Green and blue lines indicate cross-sectional cuts in (a) 
and (c), respectively. Red lines mark the estimated graphene boundaries, (c) Quantum dot formation 
at S = 0, illustrated in a cross-sectional cut of energy vs. position. Ec and Ey mark the edges of the 
conductance and valence bands. Tunnel barriers are formed by inducing a bandgap with an external 
E field while fixing Vt and V5 at a ratio that places the Fermi energy Ef within the gap. Uncompen- 
sated back gate voltage in the non top-gated regions enables charge accumulation in the dot and leads. 

Figure 2| Coulomb blockade at B = 0. (a) Conductance map (units of e^/h) of Coulomb 
blockade oscillations as a function of back gate voltage (V5) and the voltage on top gates 1 and 2 
(^12) at T = 110 mK in a four gate dot (device D4; see Fig. lb for labeling). The voltage on top 
gates 3 and 4 is fixed at Vts4 = 9.27 V. (b) 2D fast Fourier transform of (a) reveals the periodic 
structure. A peak corresponding to an oscillation spacing of ^ llmV in V5 reflects strong periodicity, 
while the appearance of higher harmonics reveals the non-sinusoidal nature of the Coulomb blockade 
peaks when ksT <^ Ec- (c) Coulomb diamonds are shown in a plot of AG/AVtu^s a function of 
1^12 and Vdc, where G is conductance in units of e^/h and Vdc is the DC bias across the electrodes. 
The voltages V5 = —10.7V and Vt34 = 9.27 V are held constant. Symmetric Coulomb diamonds 
suggests equal tunnel coupling to source and drain leads. The dot charging energy is Ec ^ 0.4 meV. 
(d) Conductance map (units of e^/h) of Coulomb blockade oscillations as a function of Vti2 and Vt34 
at fixed back gate voltage = —10.7V. 

Figure 3| Coulomb blockade at S = 5.2 T. (a) Quantum dot formation at S > 0, illustrated 
in a cross-sectional cut of energy vs. position, following the red line in Fig. lb. Tunnel barriers are 
formed using the exchange-enhanced z/ = gap, where high resistances due to quantum Hall fer- 
romagnetism are ideal for confinement. An isolated puddle of charge is created by fixing the Fermi 
energy in the top-gated regions at the middle of the 1/ = gap while allowing occupation of higher 
Landau levels elsewhere. Inset: Schematic illustration of the top gate geometry for device Dl. (b) 
Conductance map (units of e^/h) of Coulomb blockade oscillations as a function of Vti and Vt2 in a 
two top-gate dot (device D7), at fixed back gate voltage H = —15.4V and T = 110 mK. The slopes 
of the resonances indicate symmetric coupling to the two top gates, as expected for a centrally located 
dot. As top gate voltages are swept to more positive values, peak amplitude is suppressed, revealing 
moderate tunnel barrier tunability. (c) Coulomb diamonds are shown in a plot of conductance (units 
of e^/h) as a function of Vti and DC bias Vdc» where = H V and V5 = —14.4 V are fixed. 
Symmetric Coulomb diamonds suggests equal tunnel coupling to source and drain leads. The dot 
charging energy is Ec ^ OA meV (d) 2D fast Fourier transform of the boxed region in (b), revealing 
the strongly periodic nature of the oscillations and higher harmonics. 

Figure 4| Geometric control over Coulomb blockade period, (a) Coulomb blockade peak 
conductance as a function of back gate voltage observed in device Dl at B = 5.2 T. Black points 
represent data and the red line indicates a functional fit used to extract peak positions. The top gate 
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voltages are fixed at Va = 11.402 V and Vt2 = 12 V.(b) Coulomb blockade in device D2 at S = 7 
T (See Fig. Sic in the Supplementary Information for labeling). The top gate voltages are fixed at 
^1 = ^3 = 13 V and 1^2 = ^4 = 12 V. (c) Relative peak position, V{p) — Vq, plotted as a function 
of peak number p for the first 10 peaks of (a) and (b). V{p) is the position of peak p in back gate 
voltage Vb, and Vq is the position of the first peak (Vq = -15.4907 V and -17.9875 V for plots (a) 
and (b), respectively). Each black line is a plot of y{p) = Yll=o + 1) ~ ^{q)]p^ where p and 
q are peak index numbers, representing the average peak spacing for the particular data set. Inset: 
Simulated dot size versus measured size. Error bars represent the range of diameters expected for 
measured Coulomb blockade peak spacings within one standard deviation of the mean, (d) COMSOL 
simulation of density profile (in arbitrary units) for the lithographic gate pattern of device D4 for top 
gate voltage l^i = 1^2 = 12 V. Electron transport is restricted to resonant tunneling events through 
the constrictions, indicated by the arrows. 
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Figure 3 
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1. Gate defined quantum dots: device overview 

We analyze five fully suspended quantum dots with lithographic diameters, dmh = 2^ Am}^!-!^, of 
150 to 450 nm. Listed here are details on individual sample geometries and measurement conditions. 
Device Dl is a circular two gate dot with dmy, = 400 nm measured at S = 5.2 T with the valley 
polarized z/ = gap used for confinement. Device D2 is a circular dot with dm}^ =150 nm consisting 
of three main gates with a plunger gate suspended above (see Supplementary Figure Sic for label- 
ing). Measurements were conducted at 5 = and 7 T with the E field induced and z/ = gaps used 
for confinement, respectively. Device D3 is a circular four gate dot with dm}^ = 300 nm measured 
at = 7 T. Device D4 is a circular four gate dot with dm}^ = 400 nm measured at S = and 3 
T. Device D5 is a two gate elliptical dot with cross sectional lengths of 200 and 250 nm measured 
at S = 5 T. Additional scanning electron microscope (SEM) images of locally gated quantum dot 
devices are provided in Fig. SI. All measurements are conducted using standard Lockin techniques 
in a Leiden Cryogenics Model Minikelvin 126-TOF dilution refrigerator. An electron temperature of 
T = 110 mK is extracted from Coulomb blockade fits. 



2. Quantum point contacts defined by local gating 

By opening a bandgap beneath local top gates in a quantum point contact (QPC) geometry (Fig. 
S2a), electron transport is restricted to conductance through the constriction. Fig. S2b is a SEM 
image of a fully suspended QPC with 200 nm between the split gates. Beneath the split side gates, 
one may induce a bandgap by applying a perpendicular E field, meanwhile fixing the top and back 
gate voltages, Vt\ and V5, at a ratio that places the Fermi energy near the center of the gap. Application 
of voltage Vt2 to the central gate independently tunes the carrier density within the channel, enabling 
one to sweep the Fermi wavelength for a fixed constriction width. Full pinching off of the constriction 
is demonstrated in Fig. S2c. 

Figures S2d and S3 contain plots of conductance in units of e^//i versus central gate voltage (la- 
beled Vt2) in a device similar to that pictured in Fig. S2b. Measurements in a parallel field of = 5 
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T show conductance steps at every integer (e^//i spacing). The plot of differential conductance as a 
function of central gate voltage and magnetic field B shows the emergence of odd-numbered plateaus 
as spin degeneracy is broken near 4-5 T (Fig. S2d). A differential conductance map as a function of 
DC bias and central gate voltage exhibits behavior consistent with conductance quantization, though 
the features are obscured by Fabry-Perot like interference patterns (Fig. S3). 

3. Measured dot size extracted from Coulomb blockade fits 

The dot area extracted from quantized charge tunneling is given by A = ' ^H)? where 

AV5 is the back gate voltage needed to increase dot occupation by one electron and /3 is the carrier 
density induced by 1 V on the back gate. The global density to voltage conversion is obtained by 
fitting Landau level filling factors u = nh/eB from bulk quantum Hall data. The peak spacing 
AVi) = V^^^^ — V^^^^ of a given data set is extracted by fitting each Coulomb blockade oscillation to 
the conductance expression: 

where e is the electron charge, /c^ is Boltzmann's constant, T is temperature, A is peak amplitude, 
Vi)^res^ is the back gate voltage at resonance, and a = Cg/{Ci + Cr + Cg) is determined from the 
slopes of the Coulomb diamonds (Cg is capacitance to the back gate and Ci and Cr are capacitances 
across the left and right tunnel barriers, respectively) [1]. This expression is valid in the regime 
AE <C ksT e^/C, where AE is the single particle level spacing. This functional fit to the data 
is shown explicitly in Fig. 4a,b. Coulomb blockade data from an additional device is provided in 
Supplementary Figure S4. 



Device 


(3,m-^V-^ 


(AVb), mV 


Standard deviation of AVJ, 


Dot diameter, d = 2^/ A/tt 


Dl 


2.85 X lO^'^ 


14.0 


1.8 


565(+40/ - 34) nm 


D2 


3.19 X 10^^ 


34.1 


5.3 


342(+30/ - 24) nm 


D3 


3.05 X 10^^ 


68.8 


4.2 


260(±8) nm 


D4 


2.71 X 10^"^ 


10.7 


1.7 


663(+60/ - 47) nm 


D5 


2.85 X 10^^ 


13.7 


1.3 


571(+29/ - 26) nm 



Table 1 . List of measured dot sizes obtained from Coulomb blockade fits for five different sam- 
ples. /3 is the carrier density induced by 1 V on the back gate and (AVb) is the average peak spacing. 
Dot diameter is expressed as d = 2^ Aji^ for area A — • AV^). Error bars in the last column 
represent the range of diameters expected for measured Coulomb blockade peak spacings within one 
standard deviation of the mean. Results are plotted in Fig. 4a in the paper. 



4. Model 1: Cutoff defined by constriction saddle points 

We first discuss simulation of the spatial density profile for each quantum dot, as shown in Fig. 4d, 
which is relevant to the final two sections of the Supplementary Information. COMSOL Multiphysics, 
a commercial finite element analysis simulation tool, is used to model the spatial carrier density 
profile for a fixed top gate potential by solving the Poisson equation, V^V = —p/e^. First the 
lithographic top gate pattern, designed using TurboCAD, is imported into COMSOL and placed in a 
parallel plane defined 150 nm above the graphene flake. The flake is assumed to be a two dimensional 
metallic plate, which is justified by local compressibility measurements of the z/ = state yielding 
dfi/dn = 2 X 10~^^eVm^ at 2 T, which translates to a screening of 99% of the applied V5 voltage 
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by the bilayer (for dn = 1.5 x 10^ cm~'^) [2]. A fixed potential matching the experimental value at 
which Coulomb blockade oscillations appear is assigned to the top gates, while the flake is grounded. 
The perpendicular electric field component, Ez{x^y), is solved in a plane 5 nm above the flake. 
The approximate carrier density profile is given by the function f{x^y) = e • Ez{x^y)/eo up to 
a constant offset. This is because the graphene screens all in-plane electric fields. A contour plot 
representing the simulated spatial density profile in device D4 at a top gate voltage of 12 V is presented 
in Supplementary Figure S5. 

In the first modeling approach, whose results are presented in Fig. 4c, one may calculate dot size 
purely from gate geometry without relying on measured gap parameters. Assuming that charge accu- 
mulation in the quantum dot occurs when the carrier density exceeds a fixed cutoff do, the dot size is 
defined as the area bounded by the intersection with the density distribution f{x^y) with the cutoff. 
Assuming that the tunneling probability into the dot decays exponentially with barrier width, place- 
ment of the cutoff at the saddle points of the density profile within the constrictions enables maximal 
tunneling without loss of confinement. Thus, the simulated dot area computed by this method is 
simply the area bounded by the closed contour of /(x, y) at fixed density do. 

The model is used to estimate the extent to which dot size should change in response to a changing 
top gate voltage. The dot diameter, d = 2^/A/7T, for sample D2 is computed to be 170 nm at top 
gate voltage 1^ = 9 V and 173.5 nm 3iVt = 13V. This 3.5 nm increase in diameter, a 2 percent 
change, over a 4 V range is substantially smaller than the experimental error bars due to fluctuations 
in peak spacing (Fig. 4). This prediction is consistent with the overall experimental observation of an 
approximately constant dot size over the measurement ranges presented in this paper. 

5. Model 2: Cutoff determined from measured u = gap and offset from charge 

neutrality point 

In the second modeling approach, we determine quantum dot dimensions using both the simulated 
density distribution and experimental gap measurements. First the spatial density profile induced by 
the top gates, /(x, y), is modeled using COMSOL following the procedure described in the preceding 
section of the Supplementary Information. To determine the quantum dot size from the density profile, 
we impose a cutoff di determined by the width of the gap (u = or E field induced) above which 
charge accumulation begins. Explicitly, di = (y+ — V-)f3, where and V- are the positive and 
negative back gate voltages at which a plot of conductance versus V5 intersects le^ /h. Additionally 
we account for overall offsets in density due to a displacement of the measurement voltage from 
the charge neutrality point. This offset is given by doff set = {Vmeas - Vcnp)^, where Vmeas is 
the measured back gate voltage at which Coulomb blockade oscillations appear, and Vcnp is the 
measured back gate voltage at which the charge neutrality point appears at S = when the top gates 
are fixed at the potential defined in the density profile simulation. Thus, the spatial density profile with 
proper offsets included is g{x^y) = f{x^y) — dgat + doffseu where dgat is the saturating value of 
f{x^y) deep within the tunnel barriers defined by the top gates. Physically dgat is the offset in carrier 
density induced by the back gate that places the Fermi energy at the center of the bandgap. Similar 
to Model 1, the simulated dot area computed by this method is the area bounded by the contour lines 
of g{x^ y) at fixed density di. Results of Models 1 and 2 are presented in Table 2 and Supplementary 
Figure S5. 
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Device 


do, m ^ 


dl, m ^ 


Dot diameter (Model 1) 


Dot diameter (Model 2) 


Dl 


7.87 X 10^4 


6.81 X 10^4 


489 nm 


564 nm 


D2 


1.85 X 10^^ 


1.45 X 10^^ 


170 nm 


204 nm 


D3 


7.59 X 10^4 


8.34 X lO^^^ 


340 nm 


300 nm 


D4 


4.29 X lO^"^ 


3.90 X 10^4 


504 nm 


537 nm 


D5 


6.96 X 10^"^ 


6.59 X 10^"^ 


533 nm 


562 nm 



Table 2. List of simulated dot sizes obtained from Models 1 and 2 for five different samples. 
Carrier densities do and di represent the cutoff values above which charge accumulation begins in 
Models 1 and 2, respectively. Dot diameter is expressed as = 2^ A/tt for area A. Diameters from 
Model 2 were computed at the maximum value of offset doff set = (Vmeas — Vcnp)/3. Results are 
plotted in Supplementary Figure S5. 
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Supplementary figure legends 

Figure SI | Scanning electron micrographs of gate defined quantum dots in graphene (a) 

Tilted SEM image of a four gate quantum dot device, (b) SEM image of a three gate quantum dot 
taken with the Inlens detector. An additional plunger gate used to control the density in the dot is 
suspended above the lower gates. The suspended graphene bilayer is faintly visible below the gates, 
(c) Colored SEM image of a device similar to that pictured in (b) taken with the SE2 detector. Red 
lines mark the estimated graphene boundaries. The gate geometry is nearly identical to that of device 
D2, and the top gate labeling that accompanies the data in Fig. 4b of the paper is provided here. 

Figure S2| Graphene quantum point contacts defined by local gating (a) Tilted false-color 
SEM image of three suspended quantum point contact (QPC) devices in series, (b) SEM image of 
fully suspended QPC with 200 nm between split gates (labeled with blue arrows). The central gate 
(labeled with a pink arrow) tunes the carrier density in the channel, (c) Pinching off of the constriction 
at S = is illustrated. This is a plot of conductance (in units of e^/h) versus central gate voltage, 
Vt2. The split gates are fixed at Vn = -9.75 V and the back gate is at H = 11-4 V. (d) Plot of 
conductance in units of e^/h versus central gate voltage, Vt2, in a device similar to that pictured in 
part b. The split gates are fixed at V^i = — 10 V and the back gate is at V5 = 11.9 V. At = con- 
ductance steps are visible at values of 4, 6, and 8 e^/h (bottom panel), while steps energe at integer 
multiples ofe^/h in the presence of an in-plane magnetic field of 5 = 5 T (top panel). This behavior 
is suggestive of a broken valley degeneracy at B = and the gradual breaking spin degeneracy with 
increasing magnetic field. The central panel, a plot of AG/AVt2, shows the gradual emergence of 
the integer steps over a 5 T field range. 

Figure S3| Conductance through a quantum point contact at nonzero source-drain bias 

(a) Schematic of DC bias behavior for a QPC. The black lines in the bottom panel represent tran- 
sitions between conductance plateaus as a function of Vgate and Vdc- The energy diagrams in the 
upper panels show placement of the source and drain chemical potentials (fig, f^d) relative to the one- 
dimensional subbands at the locations marked with blue circles, (b) Plot of conductance in units of 
e^/h versus central gate voltage (labeled Vt2) in a device similar to that pictured in Fig. S2b. Steps 
emerge at integer multiples of e^/h in the presence of an in-plane magnetic field of S = 5 T (top 
panel). A map of AG/AVt2 as a function of DC bias and central gate voltage exhibits behavior 
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consistant with conductance quantization, though the features are obscured by Fabry-Perot Uke inter- 
ference patterns. 

Figure S4| Quantum confinement in device D3 (a) Coulomb blockade peak conductance vs. 
back gate voltage at 5 = 7 T. Top gate voltages are T^i = ^3 = H V and Vt2 = 10.5 V. (b) Peak 
spacings for the data in part a. (c) Coulomb diamonds in a plot of AG/ A V5, where G is conductance 
in units of e^/h and Vdc is the DC bias across the contacts, (d-f) Conductance(e^//i) of Coulomb 
blockade peaks vs. Va, Vt2, and Vts, respectively. Similar coupling to each top gate suggests a cen- 
trally located dot. (g) False-color scanning electron micrograph of a dot similar to D3. Voltages Vti, 
Vt2, and are applied to pins 3, 15, and 9&:10, respectively. The gates connected to pins 9 and 10 
are shorted together. 

Figure S5 1 Quantum dot size determined by simulations (a) Contour plot representing COM- 
SOL simulation of spatial density profile in device D4 at a top gate voltage of 12 V. The red line 
(indicated by the black arrow) is the contour line at the saddle points of the density profile, and the 
area bounded by the closed portion of this curve represents the quantum dot size calculated by Model 
1. Inset: Cross sectional cut of density profle along y = 0. The points of intersection with the 
cutoff (red circles) coincide with the red contour line that determines the dot area, (b) Simulated 
dot size versus measured size. Circles and squares represent areas calculated using Models 1 and 2, 
respectively. Error bars represent the range of diameters expected for measured Coulomb blockade 
peak spacings within one standard deviation of the mean, (c) Comparison between cutoffs in the two 
modeling approaches. Carrier densities do and di represent the cutoff values above which charge 
accumulation begins in Models 1 and 2, respectively. The black dashed lines are plots of 7/ = x ± 5n, 
where 5n ~ 10^^ is the density variation due to disorder in our suspended flakes [3]. 
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